Novel oscillation for the indirect coupling between magnetic nanoparticles 
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We study the prospect of using magnetic nanoparticles for the diluted magnetic seminconductors. 
The Ruderman-Kittel-Kasuya-Yosida formula is modified by explicitly taking into account the role 
of charge carriers inside the nanoparticles in addition to those in the medium. Calculations are 
done analytically. The final form of the coupling between nanoparticles is similar to the original 
formula except for additional phase terms which render a novel oscillatory feature with respect to 
the particle size, an interesting analogy to the same behavior when we vary their distance. This is 
to be contrasted to the previous approach which did not include the inner carriers and only favored 
a ferromagnetic coupling. The effect of inevitable deviations from a perfect sphere is estimated by 
the Born approximation. Our derivations can be readily generalized to finite temperatures. 



PACS numbers: 75.30Hx, 75.30Et, 75.75+a, 75.70Cn 



I. INTRODUCTION 

DMS (Diluted Magnetic Semiconductors) have been 
hailed as a potential spintronics device to integrate the 
computing power of CPU and the storage ability of the 
hard disc[l|, By intimately couphng them, we hope to 
come up with a device that can compute while record- 
ing, and vice versa. The way to prepare the DMS sample 
includes implanting transition-metal ions in a semicon- 
ducting template or using the digital doping. When the 
dosage increases, the fcrromagnetism mediated and fa- 
vored in the dilute limit by the charge carriers appears. 
Although annealing was found to increase the Curie tem- 
perature by forcing out the carriers originally trapped in 
the antisites, the optimum Tc is still only around llOK. 
There have been reports of room-temperature DMS, but 
they were often plagued by the appearance of magnetic 
clusters 0. This renders them unusable because the cou- 
pling between the electronics and magnetism we set out 
to preserve is lost. 

In this paper, we discuss the scenario^ of utilyzing 
the magnetic nanoparticles for DMS, whose movement 
is hindered by their size and is therefore less inclined to 
clustering. In the mean time, they are small enough to 
enter the superparamagnetic regime which allows them 
to carry a net moment but without a permanent magne- 
tizaion. Unlike the pestering clusters that appeared un- 
wantedly, the size and shape of these nanoparticles and 
their spacing can be controlled before we cover them by a 
semiconducting layer. Similar proposal has been studied 
in the context of Giant MagnetoresistanceQ for a metal- 
lic medium. Note that, although the underlying mecha- 
nism is expected to follow the theory by Ruderman, Kit- 
tel, Kasuya, and Yosida (RKKY), how the actual cou- 
pling between nanoparticles differs from the RKKY for- 
mula for point impurities is not clear. Previous workers 
assumed the validity of the original formula and concen- 
trated on the mathematical difficulty of integrating over 
spins in a sphere[^,0|, wire, slab, or semi-infinite planeQ. 
Discrepancies were found between their results and the 
experimental data0, [^, [13, [HI- We believe the root of 



this inconsistency lies in their failure to acknowledge the 
mediating role of the charge carriers inside the nanopar- 
ticles. This neglect is particularly fatal when the medium 
is semiconducting because, while the previous approach 
would start from a ferromagnetic coupling for all pairs 
of spins, ours encompasses the possibility of its becom- 
ing antiferromagnetic when the pair distance varies by as 
little as a quarter of 5 10 A - typical Fermi wavelength 
for metals. Most interestingly, our improvement predicts 
an oscillatory feature for this indirect coupling with re- 
spect to their size, in complement to the same behavior 
to their separation. 

When a long-range magnetic order is established 
among the nanoparticles, energy levels of all charge car- 
riers will be split by the Zceman energy. Since the in- 
ternal field is much stronger inside the nanoparticle, we 
shall neglect the splitting in the medium for the time 
being. However, it can be easily reinstated for the self- 
consistent calculation of the transition temperature for 
such based DMS. One fact that comes in handy is that 
these nanoparticles shall be sparsely spaced in the DMS, 
namely, their separation R is much greater than their 
size a. This allows us to take the asymptotic limit of 
the special functions involved in matching the incom- 
ing, reflected, and transmitted eigenfunctions across the 
boundary defined by the nanoparticle. In the previous 
approach^, 0, [13] j the effective couphng term after in- 
tegrating over all pairs of spins looks the same as the 
original RKKY formula: J{Ri2)ScS,i ■ Scs,2 with 



lim J{Ri2) = Jo 



COs(2fci?i?i2) 
(2fcFi?12)3 



(1) 



where Jq is a parameter of the system, kp is the Fermi 
momentum of the medium, and S'cff .i is the effective mag- 
netic moment for the i-th nanoparticle. The oscillation 
with respect to 2kpR remains valid. It is not expected 
that this oscillatory and decaying behavior should be to- 
tally anuUed when we include extra charge carriers in the 
nanoparticles. A naive guess is to replace the 2kpRi2 by 
2kp{Ri2 — 2a) + 2kp^n{2a) where kp^n is the Fermi mo- 
mentum of the nanoparticles. Detailed derivations later. 
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after considering the Zeeman energy, show that a more 
appropriate form is to lump the effect of these new mes- 
sengers in additional phase terms as: 



J{Ri 



Jo 



C0s(2fcFi?12 + 01 + 4>2) 
(2fcFi?12)3 



(2) 



where 0^ denotes the phase shift due to the different dis- 
persion inside the i-tli nanoparticle. 



The coefficients, Aip and Bip, can be determined by 
matching the boundrary condition as: 



Ai,{k') = {k'af 
BUk') = [k'af 



' k 

Jjjw {kaa)ni{k'a) - ji {k„a)ni+i (k'a) 
k 

ii{k„a)ii+i{k' a) - j^ji+i{k„a)ji{k'a) 



Finally, after properly normalizing the eigcnfunction in 
Eq.©, we denote it by ^k'lmp- 



II. DERIVATIONS AND RESULTS 

Due to the internal field, energy bands inside the 
nanoparticle are split by twice the Zeeman energy Vb, 
as shown schematically in Fig.([T|), and different Fermi 
momentum kpa- can be defined for each spin a. 








FIG. 1: The energy bands on the left belong to the charge car- 
riers inside the magnetic nanoparticles. They are split by the 
Zeeman energy, Vo- In contrast, only one Fermi momentum 
k'p exists in the medium on the right. 



A. Wave functions of free electrons 

Let us first calculate the indirect coupling between a 
spherical magnetic nanoparticle and a single impurity 
spin. The general form of the electron eigenfunction is 
easiest expressed in spherical coordinates with the origin 
set at its center. Symmetry allows us to separate the 
variables into 



(3) 



where the radial part consists of the spherical Bessel func- 
tions, ji and ni. Since the latter diverges at the origin, 
we have to exclude it for r < a: 



jiikav), if r < a; 

^Ai^jiik'r) + Bi„ni{k'r), if r > a 

where the three momenta are related by 

— L + 14 = Vo 

2m 2m 2m 



(4) 



B. Magnetic oscillation 

Spatial variation of the magnetization can be calcu- 
lated from the eigenfunction obtained above: 



{M{R)) ^ [*fe'/™T*fc''™T 

k'<k' Im 



. (5) 

Since the nanoparticle is assumed to be spherical, the 
azimuthal dependence drops out for M. Furthermore, 
the fact that i? ^ a in DMS allows us to Taylor expand 
and retain just the lowest order term in a/ R: 

mR))^Qik',)^^%^ + wik',t^^^^ (6) 



where Q[k'p) and W[k'p) are constants that depend on 
the parameters {k'p,a,Vo)'. 



Qik'p) = 
W{k'p) - 



B*i^Bi^ 



Ni^ Nil 
A*i^Bii+AiiBl^ A*^^Bii+AiiBl 



N, 



If 



N, 



(7) 



The normalization iV/o- = 1 / Af^ + Bf^ and Aicr, Bi^r, 
are to be evaluated at k' ~ k'p. In order to compare with 
the original RKKY formula, we rewrite Eq.® as: 



(M(i?)) «2v/Q2 + w^2. 



cos{2k'pR + (j){k'p)) 
(2fcp?p 



(8) 



where (j> = — arctan((3/VF). Figl^] demonstrates that, 
when we set the dispersion relation of the nanoparti- 
cles to be the same as in the medium and gradually 
turn off the Zeeman splitting, the energy obtained from 
Eq.® becomes indistinguishable from that of the pre- 
vious work[^, 0, This proves that the action of 
matching the boundary condition is equivalent to car- 
rying out the integration over spins inside the nanopar- 
ticle. Therefore, the coefficient, 2^02^^"^^ Eq.® 
can be thought of as an effective magnetic moment for the 
nanoparticle. Note that the phase term is due to dif- 
ferent dispersions relation for these two regions of charge 
carriers and was absent in the previous approach (gI. [tI. [13| . 

The effective moment, ^off, is plotted against the size 
of the nanoparticle in FiglJ] The oscillatory feature exists 
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FIG. 2: The coupling energy is calculated via Eq.(|8} by choos- 
ing the Zeeman splitting, Vo/E-p, to be 0.001 and 0.1 (repre- 
sented by dark and gray solid lines respectively). It can be 
seen that, as Vb decreases, the calculated value approaches 
that of the previous approach @, 0, (in dashed line) which 
neglected the carriers inside the nanoparticle but carried out 
the spherical integration over spins. The separation R is set 
at 30Af. 



for all densities of the medium carriers. This is unlike the 
previous prediction (gI. [tI. IT^. which favors the ferromag- 
netic coupling in the dilute limit according to the original 
RKKY formula. Similar oscillation is also found in the 
phase term in FigH) Expectedly, this phase is only ap- 
preciable when the particle size exceeds a quarter of the 
Fermi wavelength, Xp, to allow for a possible change of 
sign for the coupling. 
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FIG. 3: ScH oscillates with a/Xp where Xp = 2n /k'p is the 
Fermi wavelength characteristic of the medium. 



C. Correction due to a nonspherical surface 

Real samples are never a perfect sphere. We estimate 

the effect of this deviation or roughness by treating it 

as a scattering potential. In FiglSl an incoming wave 

is shown schematically to be scattered by the potential 

U due to the deviation of the boundary, r = f{9,(p), 

of a realistic nanoparticle from a perfect sphere. When 
— J- 

observed at a distance | R \ much longer than the size of 
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FIG. 4: The extra phase term, 0, due to the different disper- 
sion of the nanoparticle is shown to exhibit the same oscilla- 
tory behavior as Fig|3]when the particle size increases. 





FIG. 5: Schematic plot of a matter wave scattered by the 
irregularity, or the potential due to the deviation from a per- 
fect sphere, of a nanoparticle. This process is spin-dependent. 
Wavefunction ^j: ^{r) denotes the incoming wave, (f)cr{R) the 

outgoing wave, and R is the position of the observer. 



the potential, the scattered wave can be estimated by the 
Born approximation in the scattering theory: 



{r,R)U,{r)i^j:^{r) (9) 



where m is the effective mass of carriers in the medium, 
o-(^ incoming wave for spin a taken from Eq.(I3]), 

and (per (R) the wavefunction after the perturbation. The 
scattering potential Ua- comes from and is therefore de- 
fined by the deviation of the real rough boundary from a 
perfect sphere: 

U^ir^ = -f/i(r) = Vo ■ [e(/(0, ip)-r)- e{a - r)] (10) 

where Vq is the Zeeman energy and Q{r) denotes the unit 
step function. The unperturbed wavefunction is taken to 
be Eq.([3]), with which we can construct the Green's func- 
tion Ga{r, R). After substituting the unperturbed wave- 
function by the (j) in Eq. ^ , the effect of an irregular sur- 
face on the magnetic coupling can be estimated. For com- 
parison, the size dependence of the coupling energy sim- 
ilar to Figl5]is plotted in FiglHlfor both a perfect sphere 
and a rough surface characterized by r = O.OBAi? x sin 69. 
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It shows that not only the couphng strength receives a 
quantitative change, but the nodes of zero couphng get 
shifted. 



E/[Jo * 10' 



10 




-10 



FIG. 6: The energy of the couphng term is plotted as a 
function of the mean radius of the particle. The distance 
between the nanoparticle and the impurity spin is set at 
R — 30 Af. The dashed line is for a perfect sphere, while the 
solid line is for a nanoparticle characterised by the surface, 
r{9,ip) = 0.03Af X sin 661. 



D. Magnetic interaction between nanoparticles 

Finally, we are ready to generalize the above proce- 
dures to between two nanoparticles. Let us resort to 
physical arguments rather than matching the boundary 
condition at now two spherical surfaces, which no longer 
exhibits the rotational symmetry and therefore does not 
warrant the separation of variables. Imagine a virtual 
point spin sitting between the two nanoparticles. We 
already know how each nanoparticle interacts with the 
spin, which predicts a polarized Fermi sea like a tidal 
wave in Fig [7] according to Eq.®. Two tidal waves are 
generated by these two nanoparticles and meet to create 
interference patterns. Irrespective of where the virtual 
spin sits, we expect the best coupling to occur when the 
waves are in phase [l3|. Namely, when the difference be- 
tween {2kpRi + and (2fc^i?2 + 4'2) equals integer 
multiples of 27r. While the effective moment SeS,i and 
phase term can be estimated independently in Eq.®, 
we therefore expect the coupling energy between two 
nanoparticles to exhibit the form, E = J {R12) S cff ,i- Scs, 2 
with 



J(Rl2) — Jo 



cos{2k'pRi2 + + ^2) 

(2fc^i?12)3. 



(11) 



III. SUMMARY AND DISCUSSIONS 

In this article, we calculated the indirect coupling be- 
tween sparsely spaced magnetic nanoparticles in conduc- 
tors. Instead of quoting the original RKKY form for 
point spins, we redcrivc the formula by replacing the 
plane waves by the more suitable wavefunction which 
takes into account that the charge carriers inside the 
nanoparticles may have a different dispersion relation 
from those in the medium. We show that the matching of 
boundary condition for the wavefunction is equivalent to 
performing the spherical integration over the spins inside 
the nanoparticle except for an extra phase term. This 
phase term is crucial at rendering an oscillatory depen- 
dence for the coupling on the particle size. This is an 
interesting analogy to the same behavior on their sepa- 
ration. Our improvement is particularly relevant if the 
nanoparticles are to be utilized at manufacturing a new 
breed of DMS, because without the phase term the indi- 
rect coupling will mistakenly start with a ferromagnetic 
coupling for all spins. Although our discussions focus on 
spherical nanoparticles, formalisms are provided on how 
to deal with irregular surfaces. Our derivations follow 
the same procedures as the standard RKKY model and, 
therefore, can be as easily generalized to finite tempera- 
tures. 

For simplicity, we have used the same dispersion rela- 
tion for the carriers in the nanoparticle as those in the 
medium, except for the Zceman splitting. There are two 
comments here. Firstly, we checked the effect of a differ- 
ent effective mass for the nanoparticle in Figl51 It shows 
that, as the inner carriers get heavier, the oscillation with 
respect to the particle size becomes less periodic and with 
a shorter wavelength. Secondly, wc found that the same 
procedures in subsection IIC were also able to reproduce 
the dashed line in Fig[2l Instead of using the deviation 
from a perfect sphere, a Zeeman splitting (energy ±Vo 
for up and down spins respectively) within the radius a 
serves as the scattering potential. Apparently, there will 
be no effect whatsoever, should Vb equal zero. So the pre- 
vious results[l,0,[I3l are reproduced at the first order of 
the Born approximation, as has been done by Pogorelov 
and coworkers [1^. The phase term that is unique to our 
approach can only be obtained when we go to higher- 
order terms. 
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FIG. 8: Coupling energy between a nanoparticle and a point 
impurity is plotted as a function of the particle size. In con- 
trast to Fig[2l part of which is shown in the inset, the period 
shortens when the effective mass of the carriers in the medium 
is larger (four times than those in the nanoparticle for this 
figure). We set the distance at 7? = 20Af and the Zeeman 
splitting Vb equals one tenth of the Fermi energy. 



